Poisson equation for weak gravitational lensing 
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Using the Newman & Penrose Jl spin coefBcient (NP) formalism, we examine the full Bianchi 
identities of general relativity in the context of gravitational lensing, where the matter and space- 
time curvature are projected into a lens plane perpendicular to the line of sight. From one component 
of the Bianchi identity, we provide a rigorous, new derivation of a Poisson equation for the projected 
matter density where the source term involves second derivatives of the observed weak gravitational 
lensing shear. We also show that the other components of the Bianchi identity reveal no new 
results. Numerical integration of the Poisson equation in test cases shows an accurate mass map 
can be constructed from the combination of a ground-based, wide-field image and a Hubble Space 
Telescope image of the same system. 
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I. INTRODUCTION 



In the current age of precision cosmology, weak gravitational lensing is an important tool in understanding cluster 
mass morphology. Because all gravitating mass (baryonic or dark) influences the path of light rays on the same footing, 
^ " gravitational lensing studies are the most direct measure of mass distribution within galaxies or galaxy clusters. 

Current weak gravitational lensing studies relate the projected matter distribution to the observed gravitational 
' shearing of images through an integral relation derived in Miralda-Escude and elsewhere. Previous work by Seitz 
bi), et al. Q and others introduces a PDE approach where the gravitational shearing is related first to the gravitational 
' potential, from which the underlying matter distribution can be found. 
I A recent paper by Kling & Keith [| introduced a new PDE approach to weak gravitational lensing that directly 

^ ' relates the gravitational shear to the projected matter distribution. By examining the weak gravitational perturbations 
' on a flat background space-time in the NP spin coefficient formalism for a tetrad of constant null vectors, Kling & Keith 
showed that one component of the Bianchi identity provides a first-order, complex partial differential equation 
00 ' relating the gravitational shearing of images to the projected matter density. Using a Green's function, the authors 

, were able to show that their PDE was equivalent to the integral equation derived by Miralda-Escude 0]. 
\^ • This paper provides a more detailed derivation of the relevant PDE to first order in gravitational potential and 
using a non-constant NP null tetrad. To this order, we show that the other components of the Bianchi identity (that 
2? were not examined in Kling & Keith Q) yield no new information related to gravitational lensing. We show that a 
real Poisson equation naturally arises from this complex, first-order PDE that relates the gravitational shear to the 
projected matter density. The Poisson equation derived here is related to the Poisson equation of Seitz et al. but 
is derived from first principles of general relativity. 

We numerically integrate the Poisson equation using relaxation methods for a matter distribution representing a 
■ truncated, isothermal sphere with a core radius. These numerical integrations are used to study the feasibility of the 
method for application to observational scenarios. 



II. NP FORMALISM 



In this section, we explain the physical arrangement of an observer, sources and matter distribution that leads to 
the null tetrad at the heart of the NP formalism. We then specify Ricci and Weyl tensor components and non-zero 
spin coefficients and list the Bianchi identities in the NP formalism. 
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A. Matter distribution and tetrad 



We begin by considering an observer at the origin in a weakly perturbed, flat Robertson- Walker (RW) space-time. 
We assume that the perturbation is localized in some region far from the observer and that light sources lie beyond 
the perturbation, near the line of sight. We assume that these sources do not influence the space-time metric. For 
simplicity, we take the center of the matter distribution that perturbs the metric to lie at some distance along the —z 
axis. The space-time metric takes the form 



d~s^ = {I + 2if) dt^ -a{tY {I ^2ip){dx^ +dy^ +d7J}, (1) 

where the perturbation Lp depends on the proper distance in the background RW metric at the time light rays from 
the sources pass the lens on their way to the observer. The metric, Eq. [l]is conformal to the static metric 



ds^ = (1 + 2lp) dt^ - (1 - 2ip){dx^ + dy^ + dz^}, (2) 



and the light rays of the two metrics are identical. 

The NP spin coefficient formalism is based on a tetrad of null vectors, 



A^ = {^,n^m^m''}, (3) 

where m° and are complex spatial vectors spanning the cross section of a pencil of light rays and £° and are 
real null vectors. The tetrad vectors satisfy t°'na — 1 and m°-fha = — 1 with all other products equal to zero. The 
first null vector, is chosen to be tangent to the past-directed null geodesies. 

For our purposes, we will be considering null geodesies that connect our observer at the origin with sources lying 
beyond the lens near the —z axis. Since we are pursuing weak gravitational lensing to first order in the gravitational 
potential, we will assume that the null geodesies connecting the observer and source are those of the background 
metric 

ds^ dt^ -{dx^ + dy^ + dz^}. (4) 

Variations in these null geodesies would be first order corrections in the gravitational potential. Using the complex 
variable 

C = cot e/2 = + , (5) 

^/x + + z"^ — z 

we can write the cartesian coordinate components of null tetrad vectors as 

r = ^=_l__(-i-cCC + C,*(C-C),-i + CC), (6) 
(-i-CC,-(C + C),*(C-C),i-CC), (7) 

(0,l-C2,-z(c2 + l),2C). (8) 
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All of these vectors are null in the background metric, Eq. 21 When C = 0, the spatial part of is aligned with the 
— z axis. n° and to° are parallel propagated in the background metric along the £° congruence. Geometrically, m'^ 
is a complex spatial vector perpendicular to the spatial part of that lies in the cross section of the light bundle 
surrounding £"". 
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B. Spin coefficients and curvature tensors 



The Ricci rotation coefficients are defined by 



l^,k^\]K^aKb. (9) 

and in the NP formalism, these are arranged into twelve complex NP spin coefficients. In our examination of the 
Bianchi identity, we will assume that ^ is small so that the light rays lie close to the —z axis. The spin coefficients 
appear in the Bianchi identity multiplied by a Weyl or Ricci tensor component which we compute to first order in the 
perturbation, Hence, we will keep only the zeroth order in C, contribution to the spin coefficients. In the appendix, 
we show that the only non-zero spin coefficients to this order are 



where = + ij^ + . 

The Ricci and Weyl tensors are decomposed into tetrad components by contracting the coordinate components 
with the tetrad vectors. We assume that the gravitational perturbation and C are small, and discard terms of the 
form C, x Lp. In the appendix, we find that the non-zero NP Ricci curvature components are 



^•oo = *22 = 2$ii = (11) 



The first order Weyl tensor components are 



*4 = 






*3 = 


Ij^iVxz - 


iVyz) 


5-2 = 







(12) 

The four directional derivatives associated with the null tetrad play an important role in the NP formalism and are 
given special names: 



D=tda, A^n^da, S^m^da, 'S^fh'^da. (13) 
Since the potential is time independent, to zeroth order in Q, these derivative operators act as 



C. Bianchi identity 



The Bianchi identity can be expressed in terms of the Ricci and Weyl tensor components, spin coefhcients and the 
directional derivatives in eight complex equations 5] . A central purpose of this paper is to examine all of the Bianchi 
identity equations under the context of weak gravitational lensing. Here, we list the Bianchi identity components that 
relate the Ricci and Weyl curvature tensor components, keeping only terms which are first order ip and discarding 
terms of order (^ip. 



(15) 
(16) 
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3(^*1 - D*2) + 2D$ii - A$oo = -9p^2 (17) 

3(A*i - 5^2) - 2(5$ii = -e^i^-i (18) 

3(^^'2 - L'^'a) + 2^$ii = -Gp^-a (19) 

3(A*2 - •J^'s) + £)$22 - 2A$ii = -9Ai4'2 (20) 

_ D^4 = -p^4 (21) 

A*3 - ^^'4 + (5$22 = -4^*3 (22) 



These equations are trivial identities if one uses the definitions of the Ricci and Weyl curvature tensor components 
from the right hand sides of Eqs. [Tl] and[T2l 

Our purpose is to examine these equations as field equations, where we assume that we know the Weyl tensor 
components and wish to solve for the Ricci tensor components. To this end, we will use the equivalences of the left 
hand sides of Eqs. [Tl] andfT2l 

We note that this is the opposite of the program of Newman and Penrose, who studied gravitational radiation using 
these equations. There, vacuum space-time was studied so that the Ricci tensor was zero, and the Bianchi identities 
provided field equations to examine Weyl tensor components at infinity, which represented gravitational radiation. 



III. WEAK GRAVITATIONAL LENSING 



In gravitational lensing, the three-dimensional matter density that leads to the deflection of light is projected into a 
two-dimensional surface matter density in the "lens plane." The null geodesies connecting the observer and a source 
are replaced by two straight (asymptotic) rays that form a deflection angle at the lens plane. This leads to a lens 
equation mapping image locations in the lens plane into source positions in the source plane. The lens equation map 
depends on first derivatives (in the lens plane) of the projected gravitational potential. The Jacobian of this mapping 
involves second derivatives of the projected gravitational potential. This Jacobian can be used to map the shape of 
extended sources in the source plane into extended images in the lens plane, as is demonstrated in 1 and 0. 

We define ip{x, y) as the projected gravitational potential. 



V'(a;,y) = 




where the integral is taken from the observer position to the source position along the path connecting the source and 
observer. The projected gravitational potential depends only on the position the ray connecting the observer to the 
source strikes the lens plane (whose coordinates are x and y). 

The goal of weak gravitational lensing is to determine the projected matter density, k, from observed weak grav- 
itational "shears:" 71 and 72. These three quantities arise in the Jacobian of the lens mapping and are defined 
as 



H=\{i^xx+^yy) = ^AV (24) 

71 = 'l^{'<Pxx - i^yy) & 72 = i^xy (25) 

(Note that this matter density k is unrelated to the spin-coefficient k of Eq. IA.5l ) Using a Fourier transform technique, 
one can show that 



where 77 is the angle between the two-dimensional position vectors r and r* in the lens plane. This integral equation 
is the basis for observational weak gravitational lensing. 



IV. PROJECTION OF CURVATURE TENSORS 



The connection of our NP formalism to gravitational lensing is accomplished by projecting each of our Bianchi 
identities and the relevant curvature components along the paths connecting the observer and source. While the 
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Bianchi identity relates curvature components at each point along the light ray connecting the source and observer, 
the projection represents an accumulation similar to that studied in gravitational lensing. 

We begin by projecting the and $00 curvature components. Since we are interested in weak lensing (where the 
deflection is minimal), we project by integrating along the spatial part of the past light rays given by i"^ in the tetrad, 
Eq. [6] Since we want to work to lowest order in the potential, this is equivalent to integrating along the z axis from 
the observer to the source. Using a pre-script "L" to denote a quantity projected into the lens plane, we define 



r= 1 

L$oo = J $00 = 2 '^^^^ + ^yy^ ' ''^ 

where we use Eq. [TT] with Eq. [221 and the assumption that cp^ is zero far from the lens. Likewise, the Weyl tensor 
component ^0 is projected as 



L^-o = J dz^a^ - (V'^^ - ipyy - 2i'ip^y) . (28) 

We see an immediate connection between the curvature tensor components and gravitational lensing: l^oo = ^ and 
L^o = 7i — *72- In Kling & Keith 0], l$oo and l'^q are shown to obey an analogous integral relation to Eq. 
The remaining Ricci tensor components project as 



L*22 = 2l*ii = L*00 = K- (29) 

Likewise, = l^o = 7i + *72, but the projections of ^'i and are both zero because both involve a z derivative 
and we assume and (py are zero far from the lens. The projection of is 



-^2=1 rfz^ (</?.z-^VVj =-^(V'x. + V'ya) = -^K- (30) 



That the Weyl tensor component is equal to the mass density is consistent with the role ^'2 plays in the Bondi 
definition of the mass of an asymptotically flat space-time. 



V. PROJECTION OF BIANCHI IDENTITIES 



Each Bianchi identity, Eqs. [TSII^ can be projected in a similar manner. While each identity holds under the 
projection when using the right hand sides of Eqs. [TT] and I12| our interest is in forming field equations for the 
components of the Ricci tensor (or matter density) given the Weyl tensor (or gravitational shear) as a source under 
the conditions of gravitational lensing. This justifies throwing away boundary value terms or terms that are "small" 
under the projection relative to other quantities in the equation. 

Under the projection, 5 and 5 directional derivatives lie in the lens plane (to zeroth order in C,) and the integration 
can be passed through the derivative. This implies, for example, that the projection of S'^i will be zero since the 
projection of ^'i is zero. When the projection is applied to quantities on which the D and A derivative operators act, 
one has 



dzD^ii^-^j^ d.2^$ii«0, (31) 



where we will assume that all the curvature tensors are approximately zero far from the lens. 

For four of the Bianchi identities, Eqs. [TBI [TTl [201 and[2Tl the entire left hand side of the Bianchi identity "vanishes" 
under projection. This means that these equations have no content as field equations for projected Ricci or Weyl 
tensor components. The equality of the identity under the projection holds exactly because the small boundary value 
terms that we choose to neglect cancel the small projected term on the right hand side of the equation. 

Using Eqs.[TTl[T2land[T4l one can see that the four remaining Bianchi identities represent two independent equations; 
Eq. [THlis equivalent to Eq. [THl and Eq. [THlis equivalent to Eq. [2H Under projection, these equations become 
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<^(3l*2 + 2l$ii) 6 / dzp-^i, (32) 

J Zo 

Sl^oo-Sl'^o = 4 / ' dzp^i. (33) 



Recall that our physical setup was to have a lens at z = z; and the observer at the origin. At the observer, p diverges 
because the past light-cone which generates the null tetrad has its apex at the observer. For physically reasonable 
gravitational potentials of the form ip(\f— n|) far from the lens, the combination /5\E'i is non-zero but finite and small 
at the observer's location. In addition, because the physical gravitational potential will only be substantially different 
from zero near the lens, the combination p'i'i is a nearly odd function in z — zi due to the z derivative in the definition 
of 5*1. Thus, the projection of p^'i is very small compared with individual terms on the left hand side of Eqs.[32land 
[53] for physically reasonable models. 

Neglecting these small terms, we see that Eq. [32] ceases to be a differential equation: 



<5(3l*2 + 2l$ii) = 0^3l*2 + 2l*ii = c (34) 

for some constant c. Using the definitions of the projection, Eqs. [251 and [501 we see that the constant is zero, and this 
equation has no implications for gravitational lensing. 
The remaining Bianchi identity becomes 



S L^oo = S L'i'o, (35) 

which is a field equation for the projected mass density, l^qq or k, with the observed gravitational shear, — 
7i — il2, as a source. Eg. 1551 was the field equation reported in Khng & Keith Q, so that the different choice of tetrad 
between the two papers was not important. In Kling & Keith it is shown that Eq. [35]is equivalent to the integral 
equation, Eq. 1261 by use of a Green's function on the sphere. 

By applying the 5 operator to both sides of Eq. [35l a Poisson equation is obtained: 



This partial different equation provides a straight-forward way to determine gravitational lensing mass maps from 
observed weak lensing shear. Equation [SS] is fully accurate and consistent to first order in the gravitational potential. 
The source term, 6'^ is manifestly real when written in terms of the gravitational lensing observables 71 and 72. 
Using L^o — 7i ^ *72, we have 



5'^ L*0 = 2 (^1^^ ~ ^^yy + ^■^2a;y) - 2 il^xx - I2yy - '^llxy) ■ (37) 

Using the definitions of 71 and 72 (Eq. [25)) . one can show that the imaginary part of Eq. [571 will always be zero. 
Hence, the component of the projected Bianchi identity, Eq. 1351 informs us that projected mass density is related to 
the weak lensing shears by 



L$oo = A^K = + Kyy = 2(5^ = 713;^ - jiyy + 2j2xy (38) 

Hereafter, we will refer to the projected Ricci and Weyl tensor components l$oo and l^'o as the projected matter 
density and gravitational shears, k and 71 — ij2- 



VI. NUMERICAL METHODS 



We examine the feasibility of integrating Eq. [3H] by introducing a sample matter distribution and computing k di- 
rectly from Eq.[25]and indirectly from numerical integration. Our goal is to evaluate, under the simplest observational 
scenarios possible, whether numerical integration of Eq. [35] is a promising approach to determining accurate mass 
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maps. Keeping with the simple approach, we choose to use relaxation methods for our numerical integration based 
on Press et al. [8|]. 

In relaxation methods, one attempts to set a grid spacing that is appropriate to the rate of change of the functions 
being integrated. However, we choose to constrain the grid spacing based on approximate number counts of objects 
in the background of a lens that are resolved by the telescope in question. Observational weak gravitational lensing is 
an inherently statistical study since the natural shape of unlensed background objects is not known and only assumed 
to be on average circular. For the purposes of this paper, we assume that there are not intrinsic background galaxy 
alignments. 

Far from the center of a matter distribution, the gravitational lensing shear signal is weak, and a relatively large 
number of objects is required to form a good estimate of the underlying shear. Closer to the center of the distribution, 
a smaller number of objects is needed because the shear is substantially larger. Generally, 30-40 objects far from the 
center of the matter distribution and 10 objects in the central region are sufficient for a reasonable statistical weak 
lensing analysis. 

For the purposes of our numerical tests, we consider that one can measure average values of 71 and 72 in each box 
and assign these to grid values. While we could compute analytic derivatives of 71 and 72 to use on the right hand 
side of Eq. 1381 we will use numerical derivatives based on the grid values. Hence, as the grid spacing grows, precision 
will be lost in the source term of our Poisson equation. This is appropriate because we are attempting to model actual 
practice where the data must be binned to overcome the statistical nature of the measurements. 

For our model, we consider a truncated isothermal sphere matter distribution with a core radius, whose three 
dimensional matter density is given by 



27r G r2 + r? + ' 



(39) 



where ay is the velocity dispersion and rc and rt are the core and tidal radii of the distribution. The model is inspired 
by the truncated NFW model proposed by Baltz et al [§], but yields simpler analytic expressions for the projected 
matter density and gravitational shear expressions than either NFW or truncated NFW models. For r rt, the 
model represents an isothermal sphere with a core radius. The total mass of the model (extended to infinity) is finite 
and is given by 



Mt„t = 



G{r, + rt)' 



(40) 



and the three-dimensional gravitational potential approaches a point mass potential for large r. 
The projected mass density is 




for the radial parameter s — \/ x'^ + in the lens plane. The projected gravitational potential, ip found by integrating 
the two-dimensional Poisson equation and using the constants of integration to make the potential finite at the origin 
is 



zP 



2TT<jlr^ 



\J + r1 — \/ + rj? + Tc In 



- rt In 



rt 



+ + r1 



rt 



(42) 



The Ricci tensor or k, can be explicitly calculated from Eas.l29las 



^2 ^2 



(43) 



The gravitational shears for this matter distribution are 
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72 = 



71 = 



(44) 



for 



S 



1 



1 



(45) 



+r2(rf + v^s^ +r2)2 



_|_ ^2 (^^^ ^g2 ^ ^2)2 ' 



and cartesian coordinates {x,y} in the lens plane. We will compute 71 and 72 from Eq. 1441 and assign these values 
to each grid point. We emphasize that we will take numerical derivatives of these values to simulate the case of 
observational data. We note that Poisson's equation is independent of a^, so that the relative errors we cite in this 
paper are invariant under scaling ct„, which controls the mass of the cluster. 

For our numerical modeling, we choose to set the core radius to 250 kpc and the tidal radius to 1.5 Mpc. The total 
mass of the cluster for a velocity dispersion of 1500 km s^^ is equal to 2.1 x 10^^ solar masses. At a radius of 0.5 Mpc 
from the center, the ratio of the densities p{r = O.dMpc) / p{r — 0) is only 0.18, but 6.0 Mpc from the center, this 
ratio has fallen to 1.0 x 10~^. 



In this section, we are interested in testing the integration of the Poisson equation, Eq. 1381 under conditions 
representative of typical weak lensing measurements. We consider ground based observations by wide field telescopes, 
space based observations made by HST, and a hybrid, two-grid combination where one has both a quality ground and 
space based image. 

As our primary comparison, we look at the scaled error in the projected mass density in each box, found by 
subtracting the value of k found by numerical integration of Eq. [35] to the exact expression in Eq. 33] and dividing by 
the correct value of k at the center. (The more common relative error is less meaningful, as it approaches one at the 
edges of the grid because we set the edge values to zero as the boundary condition, even though this is in fact a small 
error.) We also compare the total integrated masses, which we estimate as the sum of the mass densities in each box 
times the box area. 

As an example of a standard ground-based, wide-field telescope usage, we consider the number counts and field of 
view (fov) for the DLS project reported in Wittman et al [l^. Here, 4-m class telescopes with 35 arc sec fov were used 
to observe five well separated fields for weak lensing signals. Filtering on well observed objects using the standard 
software for source extraction (SExtractor Tl]) and shape measurement (ellipto [i3l)j yielded 20 objects per square 
minute. Assuming that a lens is at a redshift of 0.45 with minimal dithering of the stacked images, this translates 
into an observed lens plane that is approximately 12 Mpc across in the standard cosmology. 

Hence, while one would like a smaller grid spacing for resolution, dividing the field of view into 25 x 25 boxes that 
are 0.48 Mpc across yields an appropriate average of 40 objects per box. Adding one out of field of view set of data 
points to which one can assign a boundary value of zero yields a 27 x 27 grid for relaxation of Eq. [38] We note that 
when computing numerical derivatives along the edges of the interior 25 x 25 grid, it is appropriate to use only the 
interior (measured data) points and not the added exterior edge. 

Figure [T] shows the scaled error in the projected mass density computed between numerical integration of Eq. |3S] 
and the exact value in Eq. [43] when the boundary of the 27 x 27 grid is set to zero. The error in the central bin (not 
shown in Figure[T]) is -t-31%, but the errors in other bins are all small. Along the x axis passing through the center of 
the distribution, the numerical integration slightly overestimates the mass density for values near the center, whereas 
along the diagonal, the numerical integration always slightly underestimates the correct value. When the boundary 
values are set to the correct values, which one would never know in practice, there is no significant change in the 
error estimates, implying that assigning zero to the boundary is not a significant source of error. The total integrated 
mass, estimated as the sum of the mass density times the box area, from the numerical integration underestimates 
the total integrated mass by only 8.8%. This is in part because while the central mass density value is significantly 
underestimated, the numerical integration overestimates other boxes. 

By setting the grid spacing to an arbitrarily small value, one can achieve good resolution. Using a 101 x 101 boxes 
for this field of view, corresponds to a grid spacing of 0.118 Mpc. For this spacing, the maximum relative error is 
in the central bin and is only 3.3%, and the total integrated mass is correct to 7.8%. Plots of the scaled error for 



VII. NUMERICAL TESTS 
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this number of boxes has the same general features as Fig. [T] Setting the added boundary values to the (in principle 
unknown) correct values results in extremely accurate numerical integration at this number of boxes, as one would 
expect. 

Therefore, we see that the overall limitation is the requirement of large boxes. For the scenario we have described, 
using 101 X 101 boxes results in only 2.4 objects per box. We conclude that for ground-based, wide-field telescopes 
under normal operating conditions, numerical integration of Eq.[38]by relaxation with large boxes will be good enough 
to detect weak lensing and get a sense of the overall mass, but will be insufhcient for detailed mass maps. 

We also see that the mass sheet degeneracy is not a major issue for wide field images. Setting the edges to zero 
(as opposed to the unknown correct values) does not significantly effect the results. This is because the ratio of mass 
density at the edge of a wide field compared to the center is small. 

For a space-based optical telescope such as the HST, a significantly higher number of background objects can be 
resolved, but the field of view is substantially smaller. Based on the HUDF number counts presented in Beckwith et 
al [l^, one can expect 150 objects per square arc minute for the same cuts as are used for the DLS sample. The HST 
field of view is approximately 200 arc sec across, so that for a lens at a redshift of 0.45, the lens plane is only 1.14 Mpc 
by 1.14 Mpc. However, HST can resolve approximately 1100 background objects in this field of view. Dividing the 
HST image into 11 x 11 = 121 boxes results in a grid spacing of 0.104 Mpc with approximately 13 objects per box. 

Because the mass density 0.5 Mpc from the center of a cluster is not close to zero, poor results are achieved when 
relaxing with zero as a boundary value. To see what one could achieve, we set the boundary values to the correct 
values and use relaxation to find the projected mass density. Figure [2] shows that the scaled error indicates that good 
modeling is possible if one knows the boundary values, with a maximum relative error of 3.0%. The scaled error 
curves have similar, but more gentle slopes to those of the wide-field telescopes. The total relative error in the mass 
is 1.3%. 

This modeling is unphysical because one can never know the true mass density to assign at the boundary of the 
HST image. However, if one had both a ground-based, wide field image and an HST image of the same cluster, then 
one could use the wide-field image to give interpolated values for the boundary of the HST image. 

To demonstrate the effectiveness of this approach, we assume that we have a wide field, ground based image and 
an HST image of the same system with the parameters above. We assume that the centers of these two images are 
the same and that there is no rotation between the images. After relaxing the wide field image using 25 x 25 boxes 
as before, we use the numerically obtained values of hl from the large grid to find values for the boundary of the HST 
image by interpolation. We then relax the small grid with these interpolated boundary values. 

The relative error in the k boundary values is shown in Fig. [3] along one line of the boundary of the small field 
of view grid. The points at large x correspond to the corners of the grid, and by design, the errors are symmetrical 
along each of the four edges of the small grid. We see that the approximate boundary values have a moderate relative 
error, where the relaxed values from the large grid slightly underestimate the true values near the corners of the grid 
by about 5% and substantially overestimate the values near the center of the grid by about 25%. 

Even though the relative errors in the small grid boundary values are moderate, the actual errors are small compared 
to assigning zero on the boundary of the small grid. This results in successful relaxation of the small grid in reproducing 
the mass density. The error in the central bin is less than 1% and the total relative error in the predicted mass on 
the small grid is —3.1%. (The relaxed total mass is more than the actual value.) 

Figure [4] shows the scaled errors along the x axis and diagonal as before. These scaled errors clearly track the error 
in the boundary conditions found by interpolation, as one would expect. Nevertheless, the errors are all acceptably 
small so that one can use the ground based image to provide boundary values and the space based image to provide 
accurate mapping of the mass density. 

VIII. DISCUSSION 

This paper includes a rigorous analysis of the Bianchi identity under the conditions of gravitational lensing leading 
to the derivation of a Poisson equation for the projected mass density in terms of derivatives of the observed lensing 
shear. Using the simplest possible numerical integration scheme, we show that it may be possible to use the Poisson 
equation to determine accurate mass maps when a ground-based, wide-field image and a high-resolution, space-based 
image of the same system are used in a two grid method. 

Compared with an initial investigation of the Bianchi identity 4], this paper uses a more physical tetrad that 
makes manifest the connection to the light cone of the observer. In this way, the paper is more in the spirit of the 
non-perturbative approach initiated by Frittelli & Newman [l4| and independently by Perlick [isj . The resulting 
differential relation, Eq. 1351 found in both papers is derived here in a more rigorous fashion where all the relevant 
physical approximations are made. 



10 



One might have hoped that the other components of the Bianchi identity, not examined previously in this context, 
would have yielded new differential relations between observables. In particular, one might have suspected that a 
Bianchi identity component might have contained information regarding image magnification, which can also be tied 
to the mass density. We show that this is not the case, and that no new physical information is contained in the 
remaining seven Bianchi identity components. 

A full non-perturbative treatment of image distortion is given in Frittelli et al. ,7]. It is well known that the optical 
scalars, p and a, which control the divergence and shearing of the pencil of rays connecting an extended source to the 
observer through the optical scalar equation, are related to the components of the Ricci and Weyl tensors through the 
Sach's equation. Frittelli et al. [3] show in detail how this leads to image distortion (both shearing and magnification) . 

The advance of this paper is to explicitly derive a relation between the mass density (represented by the Ricci tensor 
component l'&oo) a-nd the gravitational shears (represented by the Weyl tensor component l^o)- This connection is 
the meaningful observational relationship, but the fact that it comes directly from the Bianchi identity, rather than 
a manipulation of the optical scalar equations or the Sach's equation, would have to be viewed as a happy accident. 

The Poisson equation derived here is related to the Poisson equation reported in Seitz et al. [1], which is derived 
kinematically by considering the relationships between k, 71, and 72 and the underlying projected gravitational 
potential. When the terms that are not linear in the gravitational potential are removed from the Poisson equation 
in Seitz et al. |3|, one obtains the Poisson equation reported here. Removing the non- linear terms is consistent with 
the linearized metric used here. 

We examined reasonable observational scenarios in this paper and found that for a wide-field image alone, it would 
at least be difficult to use Eq. [35]to determine accurate mass maps. Here, we studied only relaxation methods, which 
require that the data be binned to achieve an acceptable average shear. It is possible that the more complicated to 
implement method of maximum likelihood, which can use individual data points, would result in an accurate mass 
map. 

We were able to show that the combination of a ground-based, wide-field image and an HST image could be used to 
determine an accurate mass map in the central HST region. The fact that good results were achieved using relaxation 
methods indicates that the numerical integration of Eq. [35]is a promising road for weak lensing mass maps. When 
using a wide-field image, we examined the mass sheet degeneracy and found that knowing the correct boundary values 
did not significantly improve the accuracy of the mass maps. 

Our paper did not consider several important sources of error. First, we did not add noise to the shear measurements, 
which would more accurately reflect the observational situation. The presence of small noise would make the errors in 
the two grid method larger, but we believe that our large box size (large enough to accommodate roughly 20 objects 
per box), would minimize this error. Second, we did not model missing data in our analysis. Bright nearby stars lead 
to saturation in regions on the ccd image, obscuring sometimes large regions and generally reducing the ability to 
extract background objects nearby. This would both add to the noise in the average boxed shear values and typically 
cause a handful of boxes to have no measured data. One could model these missing data points using surrounding 
data points if this was the case. Overall, we expect that these two sources of error would lead to a decrease in the 
effectiveness of the numerical integration, but would not cause the numerical integration to fail. 

Two minor issues not examined here related to the ground and HST images are not expected influence the results 
of our study. First, the two images are unlikely to be exactly centered. Second, there is typically a rotation between 
the two images. Using common stars or bright objects, one can easily orient the two coordinate systems. It is then 
simply a matter of interpolating from the large grid onto the small grid boundary. Since ground based images are so 
much larger than HST images, the center offset would have a minimal impact. A rotation would cause the relative 
error of Fig. |4]to vary in position, but the overall effectiveness should be maintained. 



IX. CONCLUSION 



This paper gives a rigorous derivation of a Poisson equation, Eq. [38l that directly relates the mass density to 
derivatives of the weak lensing shears. We show that no new information is obtained by examining the remaining 
components of the Bianchi identity, so that the full theory of weak gravitational lensing and image distortion is most 
completely described by the optical scalar equations discussed in Frittelli et al. 

Using the simplest possible integration scheme, we show that the Poisson equation can be used to detect weak 
lensing signals, but fails to give accurate mass maps for wide-field, ground based images. However, simple relaxation 
methods were shown to be successful when a ground and HST image were both present. This method of determining 
the matter distribution could prove very useful in developing highly accurate mass maps. 

Further examination of the Poisson equation is warranted. Specifically, it would be interesting to know whether 
maximum likelihood methods could yield accurate mass maps for ground based images alone. Also, the application of 
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the two grid method to mock background images that include appropriate levels of noise would solidify the usefulness 
of the method. 
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APPENDIX: NP FORMULA 

For completeness, in this appendix we briefly outline the calculation of the NP tetrad components of the Ricci and 
Weyl tensors and the spin coefficients. For the perturbed metric given by Eq. [21 the non-zero, first order Ricci and 
Weyl tensor components are given by 



i?oo = -VV Rti = -VV, (A.l) 



with V"^ ^dl + dl + dl, and 



= {^Vkk - V^<^) j ^ k Cijik = -ipjk j ^ k, (A. 2) 

where we use for time and k for spatial components, ignoring the usual index conventions. 
The NP tetrad components of the Weyl tensor are defined by 

*3 = Cafccdrn^'n^m^ *4 = -Cafccdn^m^'n^m^ (A.3) 
and the tetrad components of the Ricci tensor are 



$00 = -^Rabr£\ $10 = -ii?ab^m^ $20 = - -Rabfh'' fh'' , 



$01 - -^i^a6^m^ $11 = -i (Rabe^n'' + Rabm^'rh'') , $21 = --Rabn'^m'' 



$02 = -^Rabm^m'', $12 = -\Rabn''m\ $22 = --Rabn^n^, 

A = ^i?. (A.4) 

Given the tetrad vectors in Eqs. [B][51 one can show that the tetrad components of the Ricci tensor are independent of 
C, with $00 = $22 — 2<f>ii — (l/2)V^iy9 as the only non-zero components. The tetrad components of the Weyl tensor 
are not independent of but to zeroth order in Q are given in Eq. 1121 

The NP spin coefficients are named, complex combinations of the Ricci rotation coefficients. They are defined by 

~na-bm°-n'', tt = ~na-bm°-i'' , 
i {ia-bn'^rh'' - ma-bm'^fh'') , 

- \ {na-bi'^m'' - ma;fcm''m^) . (A.5) 



P = la-b'rrfm!', a 


= ia;bm°'m'', 


K 


/i — —na-fifh°'')Vb', A = 


-na-b'm"-m'', 


V 


e = \ {ea-bfi^i' - 


- ma;brh''&) , 


a 


l^-\{na;btn'' - 


ma-bm'^in}') , 
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We are interested in the spin coefficients for the fiat background metric. To compute the spin coefficients, it is simplest 
to change coordinates to {t, r, C, C) where the flat metric takes the form 



and the null tetrad is given by 



dr^ 



4r2 dC dC 
(1 + CO' 



(A.6) 



-1 1 



V2' V2 



0,0 



== 0,0, 



1 + CC 

V2r ' 



-1 -1 

71' 



,0,0 



to" 0,0,0, 



1 + CC 

V2r 



Although the connection for the flat metric, Eq. IA.6i is not zero, then zeros in the components of the tetrad vectors 
in the (t, r, ^, C) coordinate system make computing the spin coefficients much simpler. The non-zero spin coefficients 
are 



/3— (A.8) 
2V2r 2V2r 



Only p and fi are zeroth order in C. 
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FIG. 1: The error in the projected mass density scaled by the value of k at the origin for a wide field image with grid separation 
of 0.48 Mpc. The triangles represent the grid points along the x axis passing through the origin, while the crosses represent 
the grid points along the diagonal. 
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FIG. 2: The error in the projeeted mass density scaled by the value of k at the origin for an HST image with grid separation 
of 0.104 Mpc. The boundary of the grid has been set to the (unknown) correct values. The triangles represent the grid points 
along the x axis passing through the origin, while the crosses represent the grid points along the diagonal. 
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FIG. 3: The relative error in the projected mass density along the boundary of an HST image between the true values and 
values interpolated from a wide field image with a grid spacing of 0.48 Mpc. 
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FIG. 4: The error in the projected mass density scaled by the vahie of k at the origin for an HST image with grid separation 
of 0.48 Mpc where the boundary values are determined by interpolation from a wide field image. The triangles represent the 
grid points along the x axis passing through the origin, while the crosses represent the grid points along the diagonal. 



